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Abstract 



The phase diagram of superconducting UPt3 is explained in a Ginzburg- 
Landau theory starting from the hypothesis that the order parameter is a 
pseudo-spin singlet which transforms according to the E\ g representation of 
the D§h point group. We show how to compute the positions of the phase 
boundaries both when the applied field is in the basal plane and when it is 
along the c-axis. The experimental phase diagrams as determined by longi- 
tudinal sound velocity data can be fit using a single set of parameters. In 
particular the crossing of the upper critical field curves for the two field di- 
rections and the apparent isotropy of the phase diagram are reproduced. The 
former is a result of the magnetic properties of UPt3 and their contribution 
to the free energy in the superconducting state. The latter is a consequence 
of an approximate particle-hole symmetry. Finally we extend the theory to 
finite pressure and show that, in contrast to other models, the E\ g model 
explains the observed pressure dependence of the phase boundaries. 

PACS Nos. 74.70.Tx, 74.25.Dw, 74.20.De. 
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I. INTRODUCTION 



Currently, there is a great deal of discussion about the nature of the superconducting 
heavy- fermion compunds, especially UPt3. Much of this discussion has centered on the 
unusual nature of the superconducting state. Experiments to map out the phase diagram of 
UPt3 in the field-temperature (H — T) plane using both specific heat and longitudinal 
sound absorption |3],|| and velocity |5]|| have revealed multiple superconducting phases. In 
particular these measurements show that two superconducting phases exist even at zero 
field, as was predicted [|7| by an analysis of the free energy for a two-component order 
parameter in the presence of antiferromagnetism. |§j The resulting Ginzburg-Landau (G-L) 
theory makes additional predictions - e. g. the kink in the upper critical field when the field 
is in the basal plane. PJTD| In these theories, the order parameter transforms as one of the 



irreducible representations of the D§h point group of the crystal, either E\ or E2. [12,11 



Further evidence about the superconducting state of UPt3 comes from measurements 
of ultrasound ||i3| , [14|| and heat conduction. |15| These experiments suggest that there are 
point nodes in the superconducting gap function where the Fermi surface intersects the line 
k x = k y = and line nodes where the Fermi surface intersects the k z = or k z = tt/c planes. 
This is evidence for a d-wave E\ g order parameter which transforms like (k x k z , k y k z ). The 
theorem of Blount [ TO | states that triplet states cannot have lines of nodes when spin-orbit 
coupling is taken into account. The theorem assumes that no symmetries are present other 
than the crystal point group symmetries. It has been argued that other symmetries may 



be present in UPt3 [0 and thus lines of nodes may be present even if the Cooper pair is a 



triplet. Thus the nodal pattern may not prove singlet pairing. 

In spite of the success of the E lg model in explaining the nodal structure of the gap 
function and the existence of multiple superconducting phases, certain objections have been 
raised regarding its suitability as a description of UPt 3 . One objection is that the E lg 
theory fails to explain the isotropy of the phase diagram, or, in other words, why the phase 
diagram when the field is parallel to the c-axis of the crystal appears to be similar to the 
phase diagram when the field is perpendicular to the c-axis. [[17,18] We shall argue that 
there are similarities but also important differences and that the E\ g theory does in fact 
explain the phase diagram for both orientations of H. The other common objection to the 
Eig theory is that because it is a singlet theory it cannot explain why the upper critical field 
curve for H along the c-axis and the curve for H in the basal plane cross. ]19[ This crossing 
is maintained to be a characteristic of triplet theories alone. |20] By a careful analysis of the 
magnetic properties of UPt 3 and their contributions to the G-L free energy we will show 
that pseudo-spin singlet states can also produce this effect. 

The plan for the rest of this paper is as follows. In Sec. |TT] the overall mathematical 
approach to the phase diagram problem is discussed. It is necessary to go into the method 
in some detail: only a very careful analysis brings out the nature of the inner phase transition. 
In Sec. m we will take the free energy and use it to obtain the phase diagram when H is in 
the basal plane. The observed tetracritical point comes out in a natural way. By fitting the 
theory to the longitudinal velocity data we will obtain values for all the relevant parameters 
of our theory. Then in Sec. fV| we will obtain the phase diagram for the case when H is 
parallel to the c-axis. We will show that our theory can be fit to the data for the case 
when H is parallel to the c-axis with the same set of parameters used for the case when 
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H is in the basal plane. The near-crossing of the phase boundaries when H is along the 
c-axis is a consequence of approximate particle- hole symmetry. In Sec. we will discuss the 
magnetic properties, in particular the magnetic susceptibility. We will show the effect the 
susceptibility of UPt 3 has on the G-L free energy and how this leads to properties such as the 
crossing of the upper critical field curves for different directions of the field. In Sec. [VT| the 
phase diagram is extended to finite pressure. Finally in Sec. [VII] we make some concluding 
remarks. 



II. EFFECTIVE FIELD METHOD FOR THE PHASE DIAGRAM 

This section will be devoted to explaining the mathematical method used to obtain the 
phase diagram for UPt3 in the presence of an external magnetic field. The full problem 
is very complicated. We give first a simple example to orient the reader to the case of 
competing order parameters. The reader who is mainly interested in the overall concept, 
not the details, may read the first subsection and consult the summary figures in the other 
subsections. 



A. Simple model 

A simple system with a multicomponent order parameter and competition among the 
components is a magnet with uniaxial anisotropy. The free energy is: 

F = a 0x (T - T C )(M 2 + My) + (3 xy (M 2 + M 2 ) 2 (2.1) 
+ a 0z (T - T Z )M 2 Z + f3 z Mt + (3 XZ (M 2 + M 2 )M 2 Z . 

Suppose that T c > T z . Then at T c , the system develops a nonzero M in the x-y plane, 
its direction otherwise not determined by F. Let us say M = Mx with M(T) given by 
< M 2 >= a x0 (T c — T)/2[3 xy . The angle brackets indicate equilibrium values. The question 
we face (which adumbrates the whole theme of this paper) is: how do M y and M z behave 
below T c ? The first question has a simple answer. M y will remain zero below T c . One way 
to see this is to write an effective free energy for M y below T c by simply taking the terms in 
F which involve M y and writing the equilibrium value for M x and M z . 

F eff (M y ) = a 0x (T - T C )M 2 + 2(3 xy < M 2 > M 2 + (3 XZ < M 2 > M 2 + (3 xy M^ (2.2) 
= /3 xg < M 2 Z > M 2 + f3 xy M 4 y . (2.3) 

There will be a temperature range below T c where < M 2 >= 0. The fact that the effective 
free energy is then quartic in M y is the sign that M x and M y are degenerate, and the fact 
that the minimum of F e g is at M 2 = indicates that rotation of M in the x-y plane will 
take place only if a magnetic field (which could be infinitesimal) is applied. 
Now do the same for M z : 

F eff (M z ) = a 0z (T - T Z )M 2 Z + (3 XZ < M 2 > M 2 Z + /? 2 M 2 4 (2.4) 

= [a 0z (T - T z ) + J^a 0x (T - T C )]M 2 Z + $ Z M 4 Z (2.5) 
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There are evidently two possibilities. Either the expression in square brackets vanishes at 
positive T, in which case there is a second-order phase transition where M z appears so that 
the magnetization rotates in the x-z plane, or it vanishes at negative T, which implies that 
there is no further transition and M z = at all T. The rotational phase transition, which 
is second-order, takes place at a lower critical temperature given by 

T _ ®QzT z + ((3 xz ao x /2f3 xy )T c ^ 

An important point is that T z , the bare critical temperature for M z may be positive but T C 2 
still negative. This would be an example of the effective field suppressing a transition. If 
there is a transition (T c2 > 0), then the effective free energy is not valid for T < T c2 - it 
neglects the feedback of M z on M x . 

Of interest below will be the question of artificial terms such as 'jM^Mz in the original 
free energy. This would add a term 7[a 02 .(T c — T)/2(3 X ] 3 ^ 2 M Z to the effective free energy 
for M z . This means that M z becomes nonzero already at T c and the lower transition is 
converted to a crossover, just as if an external field in the z-direction were applied. 

For finding out whether there is a transition, whether it is second order, and computing 
the lower transition temperature, analysis of the effective free energy is all that is required. 
To find the behavior of the system below T c2 , one must minimize of the full free energy. 



B. s-wave superconductor 

We now apply the effective field method to the well-known problem of an isotropic s-wave 
superconductor to show how it works in a case which is actually quite non-trivial, but whose 
phase diagram is well understood. This system has a single complex order parameter. In 
the presence of a field, however, there are, in a certain sense, an infinite number of order 
parameters, and interesting competition among them. 

The free energy density for the system is 

/ = a (T -T c )\r ] \ 2 +P\r ] | 4 +Kj2D iV D*r)*. (2.7) 

i 

Here Di = —idi + 2eAi/hc (— e is the charge on an electron) and if we take take H in 
the ^-direction, then the gauge A = Hxy is appropriate. We have D x = —id x and D y = 
—id y + 2eHx/%c. Our problem is to minimize the free energy F — J fdV for arbitrary H 
and T. 

The method we will use is to expand the function rj(x) in a complete set of normalized 
eigenfunctions of the operator 

K(D 2 X + D 2 y ), (2.8) 

which are 

<P nk = [i/TiLlY^e-^expi-ix - k£ 2 ) 2 /2£ 2 }H n ((x - kf)/i), (2.9) 

where £ = hc/2eH, H n are the Hermite polynomials, and L y is the size of the system in the 
y-direction. We now write 
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^x) =Y J C nk <P nk {x) (2.10) 

nk 



and the free energy becomes 

F = '^J[Q.q{T — T c ) + £ n ] | C n k | + bn 1 ki,n2k2,n 3 k 3 ,n4k4C ni k 1 C n2 i i . 2 C n3 k s C n4 k 4 - 
nk niki,n2k2,n 3 k 3 ,n4ki 

(2.11) 

The coefficients in this equation are: 

AKeH 

e n = (n + l/2)^- (2.12) 

and 

b m ki,n2k2,n 3 k 3 ,n4,k4 = P J ^ x( t ) n 1 k 1 ( t ) n2k2 ( t ) n 3 k 3 ( t ) n4k4- (2-13) 

An important point is that e is independent of k and b is zero unless k\ — k 2 + k 3 — /c 4 = 0. 
We have re-expressed F as a fourth-order polynomial in an infinite number of variables C n k, 
which may be thought of formally as competing order parameters. We must minimize this 
polynomial. 

The upper critical field curve is given by noticing when the coefficient of the quadratic 
term first changes sign: 

a (T-T c )+e n = 0; (2.14) 

The highest value of H for which this equation holds corresponds to n = and the curve 

a (T - T c ) + 2KeH/hc = (2.15) 

defines the normal- superconducting phase boundary. 

Below this boundary, some but not all of the C n k are nonzero and 

C ok ~ K(T C — T) — 2KeH/hc] 1/2 = 5 1/2 . (2.16) 

This equation defines 5, which serves as a small quantity in the analysis below, the validity 
of which is thereby limited to the neighborhood of the phase boundary. 5 > in the ordered 
phase. The periodicity of the flux lattice shows that Cofc 7^ if and only if k = mq, where 



m is any integer and q = y y/3ir/£. We shall denote this condition by k G L, i. e. k belongs 
to the discrete set which constitutes the flux lattice. The discreteness reflects the fact that 
magnetic translation symmetry as well as gauge symmetry are broken in the low-T, low- 
H phase. Thus, sufficiently close to the phase boundary, only these coefficients need be 
computed and we get the familiar theory of the hexagonal flux lattice. As is well known, 
no further phase transitions take place as the field is lowered until the Meissner state takes 
over at H&. 

In UPt3, on the other hand, there is another transition when the field is reduced. Why 
does this not occur in the s-wave case? The answer is not obvious. For example, we may 
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consider the Landau level n = 1. Setting the eigenvalue equal to zero as we did for n = 
would give a critical field line with the same T c but with a slope only 1/3 of the slope of the 
H c2 curve. Why does no transition take place on this line in the H — T plane? That is, why 
is there no nonalytic behavior of the Cik on this line? What about Cok for k ^ mq ? 

To answer these questions, we must develop a picture of the effective fields present in 
the system when the symmetry has been broken. This is done by classifying different terms 
of the polynomial in Eq. |2.11 . 

Class 1: terms determing the leading behavior of Cq^, k G L. 
These are the simplest of all; the free energy is 

F = $>o(T - T c ) + e ] | C k | 2 +O(C 4 0k ) + ..., (2.17) 

k 

where only the terms relevant to the behavior of terms in class 1 have been written explicitly. 
For small 5, these terms give the simple result 

F = -5(C1) 2 + 0((C1) 4 ) (CI) ~ 5 1/2 , (2.18) 

where (CI) denote collectively the C nk which belong to class 1. For our considerations 
which are simply a matter of power counting, the indices on C are not required at this 
point. The CI are analogous to M x in the magnetic example. The conclusion is that the 
C n k are proportional to 8 1 / 2 near the phase boundary. 

Class 2: terms determining the leading behavior of C^y , k' E L'. 

We write momenta of the form k' — (m + |)g, where m is an integer, with a prime. 
Combining the Cq^ 1 builds a hexagonal lattice which interpenetrates the original one, as we 
shall see below in section [TIT] . When 5 = these variables are degenerate with the (CI) - 
these are the ones not chosen because of the breaking of the magnetic translation symmetry. 
The C2 should be compared to M y in the previous subsection. As 5 increases, they become 
less favored because they feel an effective repulsion from the (CI). The relevant terms in F 
(call them collectively Fc\,c2} are °f the form: 

Fci,C2 = £N(T - T c ) + e ] | C ok> | 2 (2.19) 

k' 

+ £ ^OAi,Ofc 2) Ofc3,O^CoA il C f o fe3 C'ofc 3 C , o fc4 + 

ki,k2,k3> k 4 

This equation repays careful examination. A first crucial point is that there are no terms 
of the form (C1) 3 (C2) or (C1)(C2) 3 . Recall that &m&i,n2fc2,ri3&3,n4&4 ^ s zero umess ki — k 2 + 
k^ — &4 = 0. However the k e L are equally spaced, so if ki,k 2 ,k^ G L, then k± G L as 
well. Similarly if ki,k 2 ,ks G L', then also k& G V . For the case k\,k 2 G L then we can 
have &3, &4 G V . Hence the only cross terms (in L and L') which survive have the form 
(C1) 2 (C2) 2 . More explicitly, 
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Fci,C2 — ^2[®o(T - T c ) + £ ] I C k' | 2 + X! ^ofei,ofe 2 ,ofe^o^C'ofciC'ofe 2 C'o^C'ofc4 (2.20) 

k' ki,k 2 ,k' 3 ,k' 4 

+ X! ^Ofci,0^,Ofe 2 ,Ofe^OA:i^Ofc2 ( ^O^C'oA:4 + CX - 
ki,k2,k' 3 ,k' 4 

+ bQk[fik' 2 fik' z fik'CQk[Ci w jjQk'j^lk' 4 + — 

y k' y k' 
In this equation 

-£>ofci,ofc 2 , 0/4,0/^ = &ofci,ofc 2 , 0/4,0/4 + &o/ci ,0/1:4,0/4,0/4 + &ofc 3 , 0/04,0/4 ,0/4 + ^ofc3,ofc 2 ,ofci,ofc^- (2-21) 

In the summations k runs over L and A;' runs over V . The idea of the effective field is to 
note that, when H < H c2 (or 5 > 0), we may write an effective free energy for the (C2): 

F eff ((C2)) = 5>„(T - T c ) + e ] I Co* I 2 (2.22) 

k' 

+ X! ^ofci,ofe 2 ,o^,o^ < C , ofc 1 Cofe 2 > Cofc^Co^ 

+ X! ^Ofci,Ofc^,Ofe 2 ,Ofc4 < CofciCofc 2 > C^/Cg^ + c - c - 
+ X! ^Ofc(,Ofc 2 ,Ofc^Ofc^Cofc(CoA; 2 Cofc^Cofc^ + ■ ■ ■ > 

y y y y 

™\ >^ 2 1^3)^4 

where the angle brackets denote equilibrium values in the ordered phase. Examination of 
this free energy is all that is required to analyze the stability of L'. Since Cok ~ S 1 ^ 2 the 
structure of this equation is 

F eff ((C2)) = -5(1 + R2)(C2) 2 + 0((C2) 4 ), (2.23) 

where R2 is a dimensionless matrix which is independent of temperature. In fact if the 
are divided into real and imaginary parts then R2 is a real, symmetric matrix. If there is to 
be no further phase transition, then all the eigenvalues of R2 must be less than or equal to -1. 
Otherwise the (C2) condense, and another lattice would form. This would mean an "inner" 
transition in ordinary s-wave materials, which does not occur. Of course the V lattice is 
degenerate with the L lattice. One may be converted into the other with application of an 
infinitesimal external current as M x may be converted to M y in the simple model by an 
infinitesimal external magnetic field. 



Class 3: terms determining the behavior of C n k for n > and k G L. 

Two subcases must be distinguished here: n = even and n = odd. If n is even then all 
possible terms come into the effective free energy for the (C3): 

F e# ((C3)) = [a (T - T c ) + (2n + I )2KeH/hc](C3) 2 + 63i < (CI) 3 > (C3) (2.24) 
+ b3 2 < (CI) 2 > (C3) 2 + 63 3 < (CI) 3 > (C2) 
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where F e ff has been written for a definite n value. The coefficient in square brackets is 
positive as long as H > H c2 / (2n + 1). 63i, b3 2 , and 63 3 are constants independent of 5 whose 
precise form need not detain us. In this field region the leading behavior is dominated by 
the term linear in the (C3): 

F eff ((C3)) ~ -(C3) 2 + 63i(Cl) 3 (C3) ~ -(C3) 2 + 63i5 3/2 (C3) => (C3) ~ 5 3/2 . (2.25) 

This resolves the question raised above. There is no phase transition at H c2 /(2n + 1) for n 
even because the coefficients determining the weight of the nth Landau level have already 
started to grow at H c2 itself. Thus the putative phase transition is converted into a crossover. 

For n = odd, only even terms in (C3) appear for parity reasons. These remain zero 
below H c2 . One argues in the same way as in case 2; since the (C3) are less stable than 
the (C2) because of the higher bare quadratic coefficient in front, the conclusion follows a 
fortiori 

Class 4 : terms determining the leading behavior of C n y, n> 0, k' G V . 

We have accumulated enough experience to write the effective free energy immediately: 
F eif ((C4)) = $>o(T - T c ) + (2n + \)2KeH/hc] \ C nk> | 2 (2.26) 

nk' 

+ 16 2J b Okl>0k2jnk t 3>n r k > 4 < C 0kl CQ k2 > C nk 'C* n , k ,^ 

ki,k2,nk' 3 ,n'k' 4 

+ 16 ^ b 0kl>ok > 3jnk2>n / k > 4 < C 0kl C 0k2 > C* nkl C* n ,y 4 + C.C. 

k\ ,A:2 ,nk' 3 ,n' k' 4 

+ bok^ok^n^y^Co^C^Cnk^Crt^ + 

By analogy with case 2, this may be written as schematically as 

F e# ((C4)) = -6(1 + i?4)(C4) 2 + 0((C4) 4 ), (2.27) 

the only difference being that R4 is a matrix in the n index as well as the k'. There are 
no terms linear in the (C4). There is no phase transition involving the (C4) - hence the 
eigenvalues of R4 are less than -1. The (C4) are always zero in equilibrium. 

Class 5: terms determing the behavior of other k values. 

It is evident from the momentum conservation condition that the effective free energy 
for k such that k (jL L, k (jL L' can contain terms such as < (CI) > (C5) 3 , for k = q/3, for 
example. There are no terms linear in the C5. The cubic terms could give rise to additional 
first order transitions in principle. It is evident that this does not occur and we do not 
consider such terms further. 

We may conclude this discussion of the s-wave case by a graphical account of how all 
transitions except the one at H c2 itself are suppressed. Fig. [l| shows the bare eigenvalue 
curves, the repulsion of some levels away from the H c2 curve and the conversion to crossover 
of others. These two effects arise from the effective field coming from the quartic term in 
the original free energy. 
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C. Effective field method for the d-wave case 



The d-wave case may be analyzed in a similar manner. For our present qualitative 
discussion we need only the fact that the order parameter becomes a complex two-component 
vector 7) = (f] x ,f] y ). The free energy again contains quadratic and quartic terms in this 
variable. It is: 

/ = a (T - T x )\ Vx \ 2 + a (T - T y )| % | 2 + ft (77 ■ rj*) 2 + ftlry • rj\ 2 (2.28) 
+ (KxDitjjD^* + K 2 D iVi D*rj* + KsD^D*^*) + K 4 £ \D zVl \ 2 

i,j=x,y i=x,y 

We have neglected certain terms which are not relevant to the present discussion. They will 
be introduced in the next section. However, we do not specialize to any particular direction 
of field, and the analysis is valid for all directions. This section generalizes the analysis 
carried out by Joynt [ETJ] for the field in the basal plane, which is the easiest case. 

The quadratic form may be diagonalized by finding the two-component vector eigenfunc- 
tions 4>nk- We then write 

•q = Y,°r*<t>nh (2-29) 

nk 



to obtain 



F = $>oCT - T cl ) + e n {H)} I C nk | 2 (2.30) 

nk 

+ ^n 1 k 1 ,n 2 k2,nzk- i ,n i k i Cn 1 k 1 C n2k2 C nA k. i C nAkA) 
nifcl ,712^2, 713^3, 714^4 

the difference with the s-wave case being that the energy levels s n (H) and the form of the b 
coefficients are far more complicated. Here T cl is the greater of T x and T y , and T c ° 2 is the lesser 
of T x and T y . The curves which are the lines a (T — T cl ) + e n (H) = are shown in Fig. ||(a). 
Crucially, however, the momentum conservation condition k\ — k 2 + k 3 — = is the same. 
We now specialize a bit to the case of UPt3. Then the solutions ao(T — T&) + e n (H) = 
fall into two classes. Half the levels have a bare T c at T ci [because e n (H = 0) = 0] and half 
at T 6 ° 2 [because e n (H = 0) = a (T cl — T c ° 2 )]. Thus when we specify the level index it must 
be stated to which class the level belongs. Call those with the higher T c (a) and those with 
the lower T c (b). Apart from this difference, the classification of states proceeds similarly to 
the s-wave case. 



1(a): terms determining the leading behavior of C na k, n = 0, k G L. 
The free energy is 

F(Co a k) = EK(T - T cl ) + e J I C 0ak | 2 +O(C 4 0ak ) + . . . , (2.31) 

k 

with the result that 

- 5(Cl) 2 a + 0((Cl)t) (Cl) a ~ (2.32) 
where 5 = -[a (T - T cl ) + e J. 
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Class 2(a): terms determining the leading behavior ofCo a k>, k' £ L' . 

Again, we write momenta of the form (m + |)g, where m is an integer, with a prime. 
These coefficients may be treated by analogy with the s-wave class 2 above. Familiar with 
the procedure, we may write down the relevant effective free energy immediately: 

F eff ((C2) a ) = -5(1 + R2 a )(C2)l + 0((C2)% (2.33) 

where again R2 a is a dimensionless matrix. This matrix is similar to R2, and we expect 
that all the eigenvalues of R2 a must be less than or equal to -1 and the usual hexagonal 
symmetry arises at H c2 . That this is actually the case has been shown by Luk'yanchuk and 
Zhitomirskii. |E2]. 

Class 3(a): terms determining the behavior of C na k for n > and k £ L. 
The analysis proceeds as in s-wave class 3. 

F eff ((C3) a ) ~ (C3) 2 a + b3 la (Clf a (C3) a ~ (C3) 2 a + b3 la 5 3 / 2 (C3) a - (C3) a ~ 5 3 ' 2 . (2.34) 

These candidate phase transitions are thus converted to crossovers by the effective field. 
Note that there is, for general field directions, no parity selection rule in the d-wave case so 
there is no distinction between n = odd and n = even. The (C3) contribute to the change 
of shape of the vortices as as the external field is reduced below H c2 , but produce no further 
phase transition. 

Class 4( a ) : terms determining the leading behavior of C na k>, n a > 0, k' G L' . 

The effective free energy is: 

F eff ((CA) a ) = -5(1 + RA a )(CA)l + 0((C4)*). (2.35) 

We again argue by analogy with s-wave case 4 that the eigenvalues of R4 a must be less than 
-1 so that the (C4) a are always zero. 



Class 5(a): other periodicities. 

Again other periodicities will not arise from states in class a, just as in s-wave, case 5. 

In discussing cases 1(a) to 5(a), we stress that we have not given explicit proofs for 
the magnitude of the eigenvalues of the different i?-matrices. It is possible to write these 
matrices formally, but they are quite complicated. However, they are very similar to the 
s-wave case, where we are certain of the result even in the absence of explicit calculation. 

Let us now turn to the the levels which start from T c ° 2 . Here we have less guidance from 
the s-wave analogy. 
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Class 1(b): terms determing the leading behavior of Co bk , k £ L. 

These terms are analogous to those in s-wave class 3. As long as H ^ 0, all possible 
terms come into the effective free energy for the (Cl)&: 

F eff ((Cl) b ) = [a (T - T cl ) + e 0b }(Cl) 2 b + b3 lab < (CI)* > (Cl) b (2.36) 
+ b3 2ab < (Cl) 2 a > (Cl) 2 b + b3 3ab < (CI)* > (Cl) 6 . 

For certain special directions of the field, the cubic-linear and linear-cubic terms may vanish, 
but we are concerned here with the general case. This leads to the result that 

F eff ((Cl) b ) ~ {Cl)l + b3 lab (Clf a (C2) b ~ (Cl)l + b3 lab 5^ 2 {Cl) b - (Cl) b ~ 5 Z ' 2 . (2.37) 

Thus these terms show crossover behavior. This fact is the apparent basis for a statement 
occasionally found in the literature that for general directions of the field there is no lower 
phase transition in d-wave systems at finite field p3|l8[. 



Class 2(b): terms determining the leading behavior of Coy , k! £ V . 
This is the crucial case so we treat it in detail. The relevant terms in F e g are: 

Fe ff ((C2) b ) = £[ad(T - T cl ) + e 0b ] \ C 0bk , | 2 (2.38) 
k> 

+ B 0aku0ak2i0bk > 3fibk > 4 < C 0akl CQ ak2 > Co^Cq^/ 

^1 

+ X! b 0aku o bk / 3t0ak2fibk > 4 < C 0akl C 0ak2 > Cq^Cq^ + c.c. 
+ b 0ak > i)0ak2)0bk3)0bk > 4 < C 0o ^C£ ofc / > C 0bk ' 3 CQ bk ^ + . . . , 

where the angle brackets denote equilibrium values in the ordered phase and 

Bo akl ,o a k2,o b k' 3> o b k' 4 = bo akl ,o ak2 ,o bk ' 3 ,o b k' 4 + &o a fei,Oi,fe4,06fe^o a fe^ (2.39) 

Since Co ak ~ S 1 ^ 2 the structure of this equation is 

F eff ((C2) b ) = [a (T - T cl ) + e 0b - 5R2 b ](C2) 2 b + 0((C2)£), (2.40) 

where again R2 b is a dimensionless matrix. The question of further phase transitions in 
this case boils down to asking whether the matrix [a (T — T cl ) + e 0b — 5R2 b ] can ever have 
negative eigenvalues. If it does, then there will be a lower phase transition. The answer 
is known for UPt3 (more formally for E representations of the hexagonal group in some 
parameter ranges) in certain limiting cases. If H = 0, then the problem reduces to a well- 
known one p|,|9|,|10|l . There are indeed two transitions and the effect of 5 R2 b is to reduce the 



bare transition temperature T° 2 (where a (T — T cl ) + e Qb (H = 0) = 0) to T c2 , the observed 
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lower transition temperature. This is the precise analog of the critical temperature for M z 
in Eq. j2.6j Also for H in the basal plane and arbitrary field strength, the problem can be 



solved [51]. At the tetracritical point, 5R2 b vanishes, so there is no effective field. Near this 
point, the effective fields can be calculated, and we will carry this calculation out in the 
next section. [cxq(T — T c {) + eo b — 8R2 b \ vanishes along a line in H-T space. This represents 
the second phase transition for this field direction. In fact all the functions involved are 
continuous and the second phase transition occurs for all directions of H for UPt3. 

Class 3(b): terms determining the behavior of C nb t for n b > and k £ L. 

The interesting new feature that arises for these terms is that we now have an effective 
field from both the (Cl) a and the (C2)&. There are more terms in the effective free energy 
for the (C3)i). However, we shall not consider this in detail, since it is evident that these 
(C3)b couple linearly to to the (Cl) a and therefore start their life at H c2 (T) where they are 
proportional to 5 3//2 . 

Class 4(b) : terms determining the leading behavior of C n k>, n > 0, k' £ L' . 

Here the relationship of terms (C4)& to (C2)& is the same as that s-wave (C3) to s-wave 
(CI). Thus there is an effective free energy of the form: 

F eS ((C4) b ) ~ [a (T - 7k) + e nb ](C4) 2 b + 0((C2)f (C4) 6 ). (2.41) 

Crossovers only are allowed for these terms. 



Class 5(b): other periodicities. 
We may neglect these for the same reasons as s-wave class 5. 

Let us now summarize the conclusions. At the Co a k condense to form a hexagonal 
lattice which we call L. At the same time, a number of other coefficients begin to grow (such 
as those in class 3(a) and 1(b)), though more slowly than the Co a k- Their growth means 
that the shape of the vortices is temperature and field-dependent, but the symmetry of the 
lattice L is unchanged. As the temperature or field is further lowered, the Co b k> become 
unstable, forming a lattice L' which interpenetrates L. This occurs by a second-order phase 
transition. This process is summarized in Fig. [| 



III. PHASE DIAGRAM: FIELD IN THE BASAL PLANE 

We have now established the mathematical method for finding the phase boundaries. In 
this section we apply the method to the quantitative construction of the phase diagram for 
the case when the field is in the basal plane of the hexagonal UPt3 crystal. We begin by 
writing down the free energy density for a hexagonal Ei or E 2 system: 
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/ = a (T - T x )\ Vx \ 2 + a (T - T,)^ 2 + ft fa ■ rf) 2 + /3 2 1 ■ rj| 2 (3.1) 
+ J] (KiD iVj D* V * + KzDirjiDffl + KzDirjjDffl) + K, £ |Z^| 2 

ij=x,y i=x,y 

+ (a eAT)(hc/2e) £ (|A%| 2 - \D i% \ 2 ) 

i=x,y 

+ a^ 2 77 • V* + a x (H 2 + H 2 )r) • 77* + a d \H ■ r)\ 2 . (3.2) 

Here 77 = (r] x ,r] y ) is the two- component order parameter, and K\, K 2 , K 3 , K4, a>o, Pi, 
P2, a x , a>z, ad and e are constants. The coupling of the staggered magnetization to 77 is 
responsible for the temperature splitting AT = T x — T y . The terms quadratic in H are 
Pauli limiting terms. They arise due to the reduction of the spin susceptibility in the singlet 
superconducting state. The effect of these terms on the phase diagram and a physical 
explanation for the relative sizes which we obtain for the various a coefficients will be given 
in Sec. |V|. The phase diagram for the case of the field in the basal plane in our theory has 
been considered previously by Joynt |21] , but we now wish to consider it in more detail and 



compare our results to experimental data. Similar calculations have also been performed in 



other models. 25.26 



It has been suggested that the spatial variation of 77 along the field direction needs to 
be considered in the calculation of the upper critical field. |1| To show that this does not 
occur, we have computed the eigenvalues for the quadratic part of Eq. |3~T] as a function of 



p , the wavevector along the field direction. The coefficient of p is positive, meaning that 
77 is uniform along the direction of the field, unless (K 2 + K^)jK\ > 3.126. As we shall see 
below, this is certainly larger than any value which can fit the upper critical field data. In 
fact the ratio is roughly unity. Let us choose a coordinate system such that H = Hx. In this 
system, the result is that we can minimize any terms in the free energy density containing 
D x by setting them to zero. Our free energy density is then 

/ = a (T - T y )|^| 2 + (P, + p 2 )\ Vy \ 4 + K,\D z ^ y \ 2 + K' 123 \D y r) y \ 2 (3.3) 
+ a (T - T x )\ Vx \ 2 + (ft + p 2 )\ Vx \ A + K 4 \D zVx \ 2 + K[\D yVx \ 2 
+ 2p 1 \r ]x \ 2 \r ]y \ 2 + p 2 (r ] : 2 r ] 2 + C.C.) 

In this equation Ki 23 = K\ + K 2 + K 3 , K[ 23 = K\ 23 — (aoeAT)(hc/2e), and K[ — K\ + 
(a eAT)(hc/2e). 

The upper critical fields of the separate components 77^ (H c2x ) and r\ y (H c2y ) are now easy 
to calculate, though the equations for the phase boundary are more conveniently expressed 
in terms of the inverse functions T c2x (H) and T c2y (H). We obtain 

T c2x = -(K + a d )/a )H 2 + H/S c2x + T x (3.4) 
T C 2y = -{a x /a )H 2 + H/S c2y + T y (3.5) 

where S c2x = -(hc/2e)(a / 'y ' K A K[) and S c2y = - (hc/2e) (a / yj K 4 K' 123 ) . S c2x and S c2y are 
the slopes of the respective H c2 curves at zero field. For a x and a d small (as we assume) 
and if T y > T x and K[ 23 > K[ then the two upper critical field curves cross. The physical 
phase boundary is the greater of the two and we obtain the well-known kink in H c2 for this 
direction of the field. We now want to find expressions for the two inner transition lines 
H* y (T) and H* X (T). 
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Consider a fixed temperature T such that H c2y (T) > H c2x (T) and ask what happens as 
the field is reduced starting from a field H > H c2y (T). As the field is lowered below H c2y (T) 
we will have rj x = but rj y 7^ 0. From the parts of the free energy involving only r] y we 
immediately see that r\ y will form a lattice. However, the lattice will be distorted from pure 
hexagonal because K4 7^ K[ 23 in general. If we choose the gauge A = —Hzy we obtain 

Vy = N y Y,c n exp{inqy - (z - nq£ 2 ) 2 / (2r 2 y £ 2 )) . (3.6) 



In this equation r 2 = y 1 K 4 / K[ 23 , N y is the (real) magnitude of r] y , and q = y \/37r(r y /£). 
Finally, c n = 1 if n is even and c n = i if n is odd. The lattice formed by |j7 y | is a centered 
rectangular lattice, athough it is perhaps more clearly imagined as a triangular lattice which 
has been "stretched" by the anisotropy. The side of the unit cell parallel to y has a length 
of 27r/g, while the side parallel to z has a length q£ 2 . 

Suppose now that we are at a temperature T such that H c2x (T) > H c2y (T) and we lower 
the field starting from H > H c2x (T). In this case as the field is lowered we will first come 
to a region where r] x ^ and r] y = 0. We may then find r) x in precisely the same way as we 
found r] y above. Using the same gauge as before we have 

Vx = N x J2 ^ exp(inqy - (z - nq£ 2 ) 2 / (2r 2 J 2 )) . (3.7) 



Here, however, we have r 2 = ^ K±/ K[ and q = y V3ir(r x /£). 

Let us return to the first case where we lower the field at constant T and H c2y (T) > 
H c2x (T) and ask what happens as the field is lowered below H c2y (T). Eventually the field 
will be low enough so that the free energy is minimized with both r] x ^ and r] y 7^ 0. 
The point where this occurs is H*(T). To calculate H*(T) formally we need to follow the 
prescription of Sec. [TT} substitute in the functional forms for r\ x and r\ y into /, compute the 
free energy F = J fdV, and determine when the coefficient of the term quadratic in N x 
changes sign. 

We must choose the functional form of r] x very carefully in this calculation for three 
reasons. First, we must allow the singularies of the rj x flux lattice to be located at different 
points than the singularities of the r\ y flux lattice while still insuring that t] x has a form 
appropriate to the gauge we have chosen. Second, we need to allow r\ x and r\ y to have 
different phases. Finally, since r\ x is arising in the effective periodic potential formed by r] y , 
we know that r] x will have the same periodicities as r\ y . (This is an application of Bloch's 
theorem HZIfl). Hence we write 



Vx = N x e ie J2 ^ exp(i(nq + z /£ 2 )(y - y ) - (z - z - nq£ 2 ) 2 / \2r 2 x £ 2 )) (3.8) 



with q = \J \/3n(r y /£). 

After minimizing the free energy with respect to the phase difference 9 the free energy is 

F = a [T - T c2y (H)\ <\ Vy \ 2 > +(P 1 + !3 2 )<\ Vy | 4 > (3.9) 
+ a [T - T c2x {H)} <\ r) x \ 2 > +(A + (3 2 )<\ Vx | 4 > 



+ 2(3, <\ Vx \ 2 \ Vy \ 2 > -2(3 2 \< vfv 2 y e 2i& >l 
< fx > + < fy > +2/hh -2/3 2 \I 2 \. (3.10) 
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Here the angle brackets denote a spatial average (< . . . >= / . . . dV), f x (f y ) is the part of the 
free energy density that depends only on rj x (r) y ), I\ =<| r\ x | 2 | r\ y | 2 >, and J 2 =< T\* x r^ j ^ %6 >. 
The inner transition (H*(T)) occurs when the term quadratic in rj x changes sign or when 

a [T - T c2x {H*)\ <\ Vx \ 2 > +Wi -2fo I h |= 0. (3.11) 
This equation contains N y which is determined by minimizing < f y > which gives us 

a [T - T c2y (H* x )] <\ r} y | 2 > +2(ft + ft) <| % | 4 >= 0. (3.12) 

We may now solve for H*(T). Once again, however, it is easier to express the result in terms 
of the inverse function T*(H). We obtain 



T:(H) = -K/a )H 2 + H/S: + T* (3.13) 



where 



(a x + a d ) - Q y a x 
l-Qv 



(3.14) 



Q* _ S c2x S c2y (l - Qy) . . 







h 






Vx 


2 >< 


Vy 


2 > 



Here 



where /3a is the Abrikosov lattice parameter: Pa =<\ % | 4> /(<l Vy | 2> ) 2 - Note that we 
may rewrite our result as 

T:(H) = T c2x (H) - [T C 2y(H) - T c2x (H)]. (3.18) 

From this equation we see that the the inner transition line for rj x is repelled away from 
the calculated upper critical field curve for r] x by an amount proportional to the separation 
between the calculated upper critical field curves for rj y and rj x . 

For the opposite case, when we consider lowering the field starting from a temperature 
T such that H c2x (T) > H c2y (T), the calculation proceeds precisely as before. The inner 
transition line, T*, for this case is related to the calculated outer transition lines as in Eq. 
3.18| so that 

T;(H) = T c2y (H) - {^Q^j [TcUH) - T c2y (H)}. (3.19) 

Here Q x is given by Eq. [3.17| . However, since rj y is becoming non-zero in the periodic 
potential formed by rj x the periodicity of the flux lattices is set by r\ x . This means that for 

this calculation q = \J y/3n(r x /£). 
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We now must calculate the integrals I± and I 2 . This is a straightforward but tedious 
exercise and we omit the details. We only note here that the results do not depend on r x 
and r y separately but only on the ratio r x jr y . This is significant for curve-fitting because 
although r x and r y cannot be obtained from the phase diagram when H is in the basal plane 
their ratio can be obtained through (r x /r y ) 2 = S c2x / S C 2y This equation follows directly from 
the definitions of the quantities involved. 

We also note that while l\ is non-zero for all values of the offset vector v = y^y + zqZ, I2 



is zero unless v is a flux lattice lattice vector or one-half of a flux lattice vector [2f| . In other 
words if tti = (2n/q)y and u 2 = (ir/q)y + q£ 2 z are the basis vectors of the flux lattice then 
I 2 is zero unless v = ^{nu-i + mu^)-, where n and m are integers. By symmetry there are only 
three distinct offset points in the Wigner-Seitz primitive unit cell of the flux lattice where 
I2 is non-zero. These are the d (v — 0), c (v = l^x), and b (v = \u 2 ) points, as shown 
in Fig. |||. Consequently, by the definition of the Q's (Eq. |3.171 ), Q x and Q y will have their 
smallest values when v is at one of these points. By the equations for the inner transition 
lines (Eqs. [3.18| and |3.19|) when Q is at its smallest the inner transition line is closest to 



the outer transition line. Hence, the inner transition line which is actually observed is the 
transition line which corresponds to the smallest value of Q. Hence, we see that the inner 
transition line must correspond to an offset vector which is at one of points d, c, or b. 

To determine which offset vector is favored and to fit our theory to the experimentally 
observed phase diagram we must calculate Q y , and from it T*(H), and Q x , and from it 
T*(H). Along with the parameters that can be obtained by requiring that the outer transi- 
tion lines fit the data (T x , S c2x , T y , S c2y , a x ) we also need P2/P1 in order to perform the 
calculation. This ratio may be determined in from the specific heat jumps at zero field at 
the outer (ACv(T y )) and inner (ACV(T* o )) transitions using 

ft _ acV(t; )/t; o , 

A &C v (T y )/T y ■ 1 • } 



From data for the specific heat jumps we obtain P2/P1 = 0.5. 

The phase diagram we obtain from our calculations is shown along with the ultrasonic 
velocity data from Ref. || in Fig. ^ along with the values of the parameters used to obtain 
it. We find that for our fit to the phase diagram the offset vector is at the b point. The fit 
is very good for the outer transition lines and H*(T) but poor for H*(T). The problems fit- 
ting H*(T) are not difficult to understand. The inner transition lines are given by equations 
such as Eq. |3.19| . These equations state that the inner transition lines are repelled from the 



continuation of the corresponding outer transition line by an amount which is proportional 
(in the case of Eq. |3.19|) to Q x /(l — Q x ). The quantity Q x has been calculated in the limit of 
very small 5, i. e. near the tetracritical point. Unlike the other quantities calculated, how- 
ever, Q x is expected to have very strong nonlinearities. The first term in Q x is proportional 
to 

,1 191 19. 

(3.21) 

This quantity is considerably less than one when the separation of the vortices is comparable 
to the core size at the tetracritical point. We find Q x = 0.333 (see below). When H < H c2 , 



< 


Vx 


2 


Vy 


2 > 


< 


Vx 


2 >< 


Vy 


2 > 
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however, the core size quickly becomes smaller than the separation and \r) x \ 2 and \r]y\ 2 are 
constant except in the region of the cores, only a small fraction of the volume. Inspection 
of Eq. |3.21| shows that Q x — > 1 and Q x /(1 — Q x ) becomes large. This magnifies the repul- 
sion between H c2y (T) and H*(T) and therefore the repulsion between the phase boundaries 
H c2x (T) and H* y (T). 

This brings in an additional error: the breakdown of the our approximation for the 
form of the order parameter. We assumed that both r] x and rj y were formed by the usual 
linear combination of lowest Landau levels. This assumption is strictly correct only at the 
tetracritical point where the inner and outer transition lines meet. It remains a reasonable 
assumption as long as the inner transition line is not repelled too far away from the outer 
transition line. These problems are not so serious for H*(T), because in this case the curve 
fits smoothly to the zero- field point, which is exact. There is no such additional constraint 
for H*(T). It is therefore necessary to incorporate some nonlinear effects in the fit to this 
line. The slope at the tetracritical point itself is corectly given by the linear calculation. 
Over the length of the line, however, we use a renormalized Q x given by fitting the slope. 
Both renormalized and unrenormalized fits are given in Fig. f|. 



IV. PHASE DIAGRAM: FIELD ALONG THE C-AXIS 

In this section we wish to take the free energy, Eq. [3.1| , and use it to compute the phase 
diagram when the field is along the c-axis (the z-direction) of the crystal. The procedure for 
finding the outer transition line (the upper critical field curve) in our theory for arbitrary 
angles of the field with the c-axis has been developed elsewhere. [^,^] We briefly review 
the procedure. 

To find the upper critical field at an arbitrary field direction one first follows the Euler- 
Lagrange prescription and demands that variations in the free energy F with respect to each 
component of the order parameter vanish. This condition gives two G-L equations which for 
purposes of finding H c2 may be linearized. The linearized G-L equations may be viewed as 
a Schrodinger equation for 77. This defines an effective hamiltonian which is a 2 x 2 matrix 
in the components of 77. One then defines a new coordinate system with one axis along the 
field and the other two axes perpendicular to it. 

It is easy to show that the component of the D operator along the field (-Di) commutes 
with the components in the other two directions. Hence Di commutes with the effective 
hamiltonian and we may rewrite any terms containing D\r\ as p\r\ where p\ is a c-number. 
When the field is in the z-direction the only terms which result from this substitution are 
terms proportional to p\, which are minimized by setting pi = Therefore in this case, 
as in the less obvious case when the field is in the basal plane, one may simplify the G-L 
equations by setting Ditj = 0. Since this procedure may be done both when the field is in 
the z-direction and in the seemingly least favorable case when the field is in the basal plane 
it is reasonable to assume that it may be done for any angle the field makes with the z-axis. 

One then defines raising and lowering operators D± = t(rD 2 ± iD^/r) /y/2 and r]± = 
(rj x ± irj y )/2. Here r is a function of the angle the field makes with the c-axis and is chosen 
to simplify the G-L equations as much as possible. One can then rewrite the G-L equations 
in terms of these quantities and expand r) + and r\- in terms of the states | n >: 
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i] + = Y, a n\n> i]_ = Y,h\n> . (4.1) 

n n 

Here D + D_ | n >= n \ n >. The states | n > are quasi-Landau levels. The problem 
then splits into finding the eigenvalues of an infinite tri-diagonal matrix. From the lowest 
eigenvalue one can then compute H c2 . 

Finding the inner transition line near to the upper critical field is also a linear problem, as 
the analysis of Sec. |H| demonstrated. To find the line rigorously we would have to calculate 
the effective free energy for all of the eigenfunctions due to the presence of the eigenfunction 
with the lowest eigenvalue, as outlined in Sec. ||. However, we have seen that the only 
transition line which is not destroyed (that is, either converted to a crossover or repelled 
to non-physical fields and temperatures) by the coupling to the lowest eigenfunction is the 
line which originates at the inner transition temperature and corrresponds to a flux lattice 
shifted from the flux lattice formed by the lowest eigenfunction. The full effective field 
matrix therefore contains levels which are pushed to unphysical fields (pushed up to high 
energy in the quantum-mechanical analogy), or have small magnitude (~ 5 3 ^ 2 ). Thus, rather 
remarkably, it will be a very good approximation to compute the inner transition line using 
only two levels. As in the case when the field is in the basal plane we then have a correction 
to the bare inner transition line which is proportional to the separation between the bare 
inner transition and the outer transition in order to find the actual inner transition line. 
Our formula for the inner transition line Ti nner (H) in terms of the bare inner transition line 
Tb are (H) and the outer transition line T 0Uter {H) is then 

Ti nn er(H) = T bare {H) — g[T 0Uter (H) — T bare (H)}. (4.2) 

When the field was in the basal plane we were able to compute the coupling constant 
g. For the field along the c-axis this computation, though straightforward in principle, is 
exceedingly complicated. Accordingly, g is found by fitting to the data. 

A key feature of the bare inner transiton line comes to light upon examining the matrix 
used to find it. This matrix is 

/ (2K 1 + K)-K' -a l 2 AT(2eH - 1) \ 

^ -a £ 2 AT(2eH - 1) (2K X + K) + K' )' 1 ' 

Here K = K 2 + K3 and K' = K 2 — K%. Note that the off-diagonal terms will vanish when 
H = l/2e. This means that if K' ph 0, as is expected from particle-hole symmetry [25|, the 



two eigenvalues will be nearly degenerate at this H and the outer and bare inner transition 
lines will nearly touch. This cancellation between the derivative terms in the free energy (Eq. 
|3.1|) proportional to K 2 and K% and the terms which couple the staggered magnetization 
(through AT) to the derivatives means we obtain an apparent tetracritical point - the two 
lines come close but do not quite touch. There are only two superconducting phases when 
the field is in the ^-direction or indeed for any direction except in the basal plane. The 
fact that K' depends sensitively on the impurity density has interesting consequences. The 
miminum separation will depend on this density. Unfortunately, the sharpness of these 
transitions also depends on the impurity density. 

The phase diagram we obtain from our calculations with the field in the ^-direction 
together with the values of the parameters used to obtain it and the ultrasonic velocity data 
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from Ref. || is shown in Fig. [5[ Note that we use the same parameters as for our fit for when 
the field is in the basal plane along with some additional parameters. As was the case when 
the field was in the basal plane the fit is very good except for the high-field, low-temperature 
part of the inner phase boundary where the linear theory is expected to break down. This 
happens because of the renormalizations discussed in the previous section. 

A virtue of the theory given here is that a striking difference between the phase diagrams 
for the two directions of the field receives an explanation. The upper critical field curve is 
smoother for field along the c-axis, and the inner transition line is much smoother. This can 
now be seen to result from the 'hybridization' of the two curves for this case, the presence 
of the off-diagonal matrix elements in Eq. [4.3| . This is absent for the other field direction, 
when the two components decouple. 



V. MAGNETIC PROPERTIES OF UPt 3 

In this section we wish to discuss the origin, effect, and relative sizes of the Pauli limiting 
terms in the free energy density. These are the terms proportional to H 2 t] 2 in Eq. |3.1| . In 
order to do this, we require some preliminary background about magnetic properties of the 
normal state. The first fact to appreciate is that the magnetic susceptibility Xij °f the 
normal state is enhanced by roughly the same factor as the mass. Because Xij is large, the 
Pauli limiting effect of the field on superconductivity is likely to be appreciable. The second 
important point is that the susceptibility is anisotropic, and the temperature dependences 
of the components are different. This is clear from the plots of the susceptibilities Xxx(T) 
and XzziT) in Fig. XxxiT) > XzziT) at all temperatures. 150] At high T, both functions 
take on the local moment form x ~ while each goes to a finite constant, characteristic 
of Pauli or van Vleck behavior, at low T. In addition, Xxx{T) has an anomaly around 15 K. 

Let us first take a theoretical approach to understanding the anisotropy in Xij- Our 
basic assumption is that UPt3 is a Fermi liquid at temperatures just above the critical 
temperature. Then the starting point is the single-particle states calculated in band theory, 
which account very well for the Fermi surface. ]32] The states near the Fermi surface are 



predominantly derived from uranium 5f orbitals with j = 5/2, as would be expected from 
Hund's rules for an actinide system with a 5f occupancy near 2. In the isolated atom, 
the j = 5/2 level is 6-fold degenerate. In the hexagonal crystal field, there is an effective 
Hamiltonian at the T point which splits the six-fold degenerate state into three doublets at 
the T point: j z = ±5/2, j z = ±3/2, and j z = ±1/2. This means that UPt 3 is likely to be an 
example of a system in which the magnetism is Van Vleck-like in the plane and Pauli-like 
along the c-axis, which is expected to be a general feature of hexagonal U-based systems. 



Let us briefly review the reasons for this expectation. If we apply a magnetic field, there 
will be both a Pauli (intraband) and a Van Vleck (interband) contribution to the suscepti- 
bility. The former is of order (g e fffi B ) 2 N(ep), while the latter is of order (g e jj{i B ) 2 /\B h \. Here 
g e ff is an effective g-factor for the coupling of the field to the total angular momentum of 
the band or bands involved. It is a dimensionless number of order unity. The Lande factor 
for £=3, s=l/2 , and j=5/2 is 6/7. Bh is the separation between the bands and N(ef) is 
the density of states at the Fermi energy. The Van Vleck susceptibility is given by 
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Xu = 2n/4 £ l<"l L ; + 2 ^>l 2 /a (l - /,). (5.1) 

Here f a , fp, E a , Ep are occupation factors and energies of the states a and (3. In view of 
the greater multiplicity of the interband transitions, we expect the Van Vleck susceptibility 
to be very important - indeed it very likely dominates the total. A band calculation which 
explicitly computes the two components reckons the Pauli contribution at 15-20%, [34] in 
rough agreement with this multiplicity argument. 

If H is along the c-axis, then the relevant matrix element (with % = 1) is: 

| < a\L z + 2S g \/3 > | 2 = (36/49)j% i/3 . (5.2) 

At the T point, states of different j z do not mix and the perturbation introduced by H is 
diagonal. The occupation factors then imply that the Van Vleck susceptibility is zero for 
this direction. If H is in the x-direction, the corresponding expression for the square of the 
matrix element is 

| < a\L x + 2S X \(3 > | 2 = (36/49)(5/2 - j z )(5/2 + j z + 1) (5.3) 

if the states a and j3 differ by one unit of j z and is zero otherwise. The Van Vleck suscep- 
tibility comes from four distinct pairs of states:(j z = —5/2, —3/2), (—3/2, —1/2), (1/2, 3/2) 
and (3/2,5/2), whenever one of the pair is occupied and the other unoccupied. The Pauli 
contribution to Xxx, on the other hand, comes only from the pair (—1/2, 1/2) when this state 
is occupied. A sheet of the Fermi surface will have an isotropic partial Pauli susceptibility 
(Xzz/Xxx ~ 1) ^ different j z values are well mixed in the wavefunction, but will be anisotropic 
otherwise: j z =l/2 implies Xzz/Xxx << 1> an d jz=3/2 or j z =5/2 implies Xzz/Xxx >> 1- As 
we shall see below, it is the anisotropy of the Pauli contribution which is critical for under- 
standing the phase diagram. This means that the central question is: what is the j z content 
of the Fermi surface, and how much mixing of different j z s is there? Band calculations give 
a clear answer to this question. They show that the parts of the Fermi surface near the 
T point and K point are predominantly j z =3/2 or 5/2, [[35|,[3lJ while the parts near the A 



point are well mixed. This is illustrated in Fig. |7|. Hence we expect a contribution to the 
Pauli susceptibility which satisfies Xzz/ Xxx ^ 1 from the parts near Y and K, representing 
roughly half the total density of states at the Fermi surface, and a contribution satisfying 
Xzz/Xxx ~ 1 f° r the rest of the Fermi surface. In treatments which go beyond band theory 
to discuss many-body renormalizations, it is found that the Pauli and Van Vleck parts are 
enhanced by similar factors | |37| |. 

Summing up these theoretical considerations, the magnetic susceptibility of UPt3 is likely 
to be dominated by interband (Van Vleck) contributions. This is particularly true for Xxx, 
which means that the anisotropy in the oberved susceptibility (x X x > Xzz) most likely stems 
from interband contributions. The Pauli susceptibility, on the other hand, is more likely to 
satisfy the opposite inequality Xxx < Xzz 

Experimentally, it is not easy to distinguish the Pauli and Van Vleck contributions to the 
susceptibility. The most straightforward way, in principle, is to measure the imaginary part 
of the susceptibility with neutron scattering. The Van Vleck contribution has a gap at low 
frequencies, while the Pauli part does not. For the present case, however, we also need to 
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distinguish the different components of the susceptibility tensor. This means that polarized 
beam experiments are required, with the associated lower counting rates. Finally, we are 
interested here in the uniform susceptibility, which means small-angle scattering. Thus this 
definitive experiment may be difficult to perform. 

A more indirect but still informative test arises from the observation that the Pauli 
susceptibility depends on the density of states at the Fermi energy whereas the Van Vleck 
susceptibility depends on a joint density of states. The Pauli part is therefore directly 
comparable to Cy/T, where Cy is the specific heat. In this regard the peak in XxxiT) at 
T=15 K |3(J (see Fig. []) is of interest. This peak is absent in the smooth curve for Xzz(T), 
and in the the specific heat Cy(T), |JT| This is consistent with the idea that the physical 
origins of Xzz and Xxx are different, and that the density of states at the Fermi level largely 
determines \ zz but not Xxx- Thus experiments, to the exent that we have them, confirm the 
theoretical picture. 

The importance of these considerations for the superconducting state is simple. 



Superconductivity affects the Pauli susceptibility in a drastic fashion. For a singlet state 
such as Ei g , the Pauli term xfj{T) is reduced to zero at zero temperature because it takes a 
finite amount of energy to break a pair and magnetize the system. Superconductivity should 
have no effect at all on the Van Vleck term, and conversely. The difference in free energies 
between the normal and superconducting states in a field is 

Fmagnetic = ~~ ^2 ^XijHiHj. (5.4) 
1 ij 

Here = xfj — X% where xfj and x% are the Pauli susceptibilities in the superconducting 
and normal states, respectively. Just below the superconducting transition we know that the 
change in the susceptibility is quadratic in 77. Hence we add to the usual superconducting 



free energy the last three terms of Eq. |3.1| which are quadratic in both rj and H. 

From the arguments above we expect that a x and ad will be smaller than a z since we 
anticipate that Xxx < Xzz- From our fits to the phase diagrams for the two directions of 
the field we find that a z is slightly more than twice a x + a^, in agreement with the physical 
picture of the susceptibility. The differences in the sizes of the a terms affect what happens 
to the upper critical field curves for the two directions of the field at high fields. At high 
fields the Pauli limiting terms in the free energy, which are proportional to H 2 , dominate 
over the rest of the free energy, which gives a contribution to H c2 proportional to H. Because 
a z > a x + a<i, H c2 when the field is along the c-axis curves down more than H c2 when the field 
is in the basal plane. Consequently, the two upper critical field curves for the two directions 
cross. This crossing is shown in Fig. ^[ We have therefore shown that the objection to 
the Ei g model on the grounds that it cannot explain the crossing of the upper critical field 
curves is invalid. 



VI. PRESSURE EFFECTS 

We have offered a comprehensive description of the phase diagram of UPt3 in the H — 
T plane. However, because of the rather large number of parameters in the Ginzburg- 
Landau free energy, this analysis is not yet sufficient to distinguish the E lg picture from 
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competing pictures such as the E 2u picture and mixed representation pictures. Consideration 
of pressure effects will allow us to do this. We will show that only E\ g is consistent with these 
experiments. The analysis in this section is an elaboration of earlier work. |39|] It is somewhat 



surprising that pressure experiments are so crucial for understanding the symmetry of the 
order parameter. Under normal circumstances, accessible pressures have only a small effect 
on superconducting parameters and qualitative conclusions are difficult to draw. In the 
present case, however, moderate pressures destroy antiferromagnetism, which restores the 
full hexagonal symmetry of the crystal structure. It is this singularly fortunate circumstance 
which makes pressure such a very powerful tool in unraveling the order parameter symmetry. 

Qualitatively, the facts are these. The magnetization disappears at a critical pressure of 
about 3 kbar. The splitting in T c also disappears at the same pressure. This shows that 
it is indeed the magnetization which splits the transition, as originally predicted M. The 
coincidence of the pressures at which these events take place rules out mixed representation 
theories such as the A-B theory ||18|| . In such theories the original splitting is due to an 
accidental degeneracy and is not related to the magnetization. 

Our aim is to understand quantitatively the phases of UPt3 in the entire (H, P, T) space. 
However, in order to understand the restoration of crystal symmetry, we first focus on the 
(H = 0,P, T) plane, so that complications due to the gradient terms can be treated sepa- 
rately. The expression for the free energy density of the coupled magnetic-superconducting 
system is then f = f s + f M + f SM , where 

f M = a M (P,T)M 2 + (3 M M* (6.1) 
f 3 = a s (P, T)ry • 77* + ftfa ■ rj*) 2 + h\l • v\ 2 (6.2) 

f SM = b\M-r 1 \ 2 + b'M 2 r 1 -r 1 * (6.3) 

We have assumed, as is conventional, that the presure dependence of fourth-order coefficients 
is weak and can be neglected. 

J'm, the magnetic part of the free energy, entirely determines the behavior of the magneti- 
zation above T c \. (Recall that T c \ is the higher of the two observed transition temperatures.) 
The experimental data from neutron scattering measurements of M 2 (proportional to the 
magnetic Bragg scattering at the (1, |,0) point) are sufficient to determine the parameters. 
At P = and T > T cl = 0.5 K, M 2 is a linear function of T/v — T, where T/v = 5 K is the 
Neel temperature. §,|tj One finds a M (P = 0, T)/f3 M = (1.6 x 1(TV|/K)(T - T N ). 

As to the P dependence, it is found that T/v is nearly independent of pressure from 
P = to P = 2 kbar and that M 2 ~ (P N - P) for T < 2 K gtj, where P N « 3 kbar 



is the critical pressure at which magnetism disappears. From the point of view of this 
paper, which concentrates on the superconducting regime T < T&, we may therefore take 
a M = a° M (P - P N )(T - T N ) « -a° M T N (P - P N ), where a° M /{3 M = 5.3 x lCr 5 /4/K-kbar. 
Note that this value and the coefficient of the expression for «m(P = 0,T)//3m have been 
corrected from an earlier paper written by one of us (Joynt) p9| . 

The pressure dependence of as(P,T) may be obtained if we assume that as(P,T) = 
a sr{T — T®) + «5pP, so that acs(P, T) is a linear function of P. ocst(T) is the zero pressure 
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and the pressure 
—11 mK/kbar in 



value of acs(P, T) which has already been determined. For P > P^, M 
dependence of T c is entirely due to the coefficient otsp- Since dT c /dP 
this region pfl|| , we find asp = o:st(11 mK/kbar). 

At P = and T < T cl , there is a competition between the purely magnetic terms and 
the coupling term fsu- Because T7 ■ 77* ~ T c \ — T for T < T cl and 77 = for T > T c i, the 
coupling term predicts that ikf should have a kink at T c \. The magnitude of the kink may 
easily be computed using Eqs. |5T2] , 
M and rj = rjx 



6~3l. Differentiation leads to two linear equations for 



r, 2 = [a 5 (T c ° - T) - (b + 6 / )M 2 ]/2(/? 1 + (3 2 ) 



(6.4) 



M 



[aA/(T w -T)-(6 + 6')^]/2/3M, 



(6.5) 



which give the behavior of the order parameters below T cl . Above T cl we have simply 



M 



o:m{T n ~ T)/2j3 M 



The slope is discontinuous at T cl : 



(b + b 



!\2 



dM 2 



dT 



dM 2 



T<T C _ 



dT 



+ 



T>T cl 



asr(b + b') 



(6.6) 



(6.7) 



If we take the approximation that the coupling (b + b') is small, then we may write the 
discontinuity as 



dM 2 a ST {b + V) 



(6.8) 



In these formulas if is assumed to be parallel to M . If these two vectors are perpendicular, 
then b should appear instead of (b + b'). The kink is observed experimentally, |§ which again 
confirms that the splitting of the superconducting transition is due to fsM- These formulas 
assume that there is only one component of M, contrary to the idea of Blount et al. |23[ 



that the moment rotates at T c \. Recent experiments have indeed ruled out the possibility 
of rotation 



We now wish to calculate the phase diagram at finite pressures, assuming that the only 
pressure dependence comes from as and olm- All other parameters are taken to have their 
zero pressure values. The only dependence on pressure in our theory of the phase diagram 
is through the quantity AT. We calculate AT at various pressures by taking T X [P = 0) and 
T y (P = 0) from our zero pressure fit, dT y /dP (recall T y > T x ), P^ and dT c /dP (for P > P/v) 
from experimental data |^T|. From T X (P = 0) and P/v we can then find dT x /dP. In Fig. |9| 
we plot the phase diagram in the H — T plane at various pressures for H in the basal plane 
using the renormalized Q x (see section PTT| ). The behavior with pressure is easily understood 
qualitatively. The main effect of pressure is to close up the splitting of the zero-field critical 
temperatures. Thus the tetracritical point moves down toward the T-axis and disappears, 
as does the A (low field, high temperature) phase. Thus the C (high field, low temperature) 
- B (low field, low temperature) phase boundary is very sensitive to pressure as it ends 
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at the tetracritical point. This is observed experimentally. [[|3] On the other hand, the N 
(normal state) - C boundary (upper part of the H c2 curve) is not very sensitive to pressure. 
Again, the agreement between theory and experiment is very satisfactory. It is difficult to 



compare these predictions with the threee-dimensional phase diagram of Boukhny et al. [44 



quantitatively. The pressure dependence of the critical temperatures given by these authors 



is not in good agreement with that of Trappmann et al. |^T[ j which we used in plotting the 
figures. The behavior of the boundaries is quite sensitive to this dependence. Nevertheless 
there appears to be very satisfactory qualitative agreement between theory and experiment, 
with one exception. The experiment shows that there is an additional phase boundary in 
the P-T plane when P > P^. This cannot be a pure superconducting transition in a two- 
component theory. We believe this to be a mixed magnetic-superconducting transition, so 
that this boundary is essentially an extension of the magnetic phase boundary. The signal 
in the sound velocity is very small. It may be larger than in the normal phase because of 
the coupling to the superconducting order parameter which is serving as a secondary order 
parameter in the transition. 

Let us compare this behavior to the behavior of the phase boundaries in the E 2u theory 



in which K 2 ~ K3 ~ 0. The best fit with this constraint is given in Fig. |10| This picture is 
in qualitative disagreement with experiment. Again, the qualitative reason for this is easily 
understood. In the E 2u theory, the difference in slope between the H c2 curves for r] x and r) y 
is due only to their differing energies in the presence of the magnetization. Once P > Pn, 
this difference is gone and the two components have identical free energies and identical 
slopes. The E 2u theory says that the N-A and A-B boundaries must move together, not 
apart, under the influence of pressure. This is in conflict with experiment. 



VII. CONCLUSION 

The Ginzburg-Landau theory is a very powerful tool in the physics of unconventional 
superconductivity. We have pushed the theory to obtain as much information as possible 
about the phase diagram. Mathematical difficulties arise when a magnetic field is applied, a 
circumstance which has made the theory of the phase diagram of UPt3 proceed more slowly 
than might have been expected. The method of classifying terms according to their behavior 
in the effective field appears to have solved the linear problem in principle, though explicit 
calculations for a general direction of the field still appear daunting. We have limited our 
treatment to the two high symmetry field directions. Most experiments are also limited to 
these directions. 

Consistent application of the method, taking into account the Pauli-limiting effect, gives 
very good agreement between theory and experiment for the E± g theory. It would be de- 
sirable, however, to have an explicit calculation of the nonlinear renormalization factors 
entering the repulsion of the phase boundaries; obtaining this by a fit, as done here, is not 
truly satisfactory from the theorist's point of view. Thr peculiar phenomenon of the H c2 
crossing is interpreted here as arising from an interplay of intraband Pauli magnetism and 
interband Van Vleck magnetism. While the picture of the anisotropic susceptibility which 
emerges is a natural one, it would be good to have some independent confirmation of it. 

The surprise of the past several years is that pressure experiments have been able to play 



24 



a critical role in sorting out the nature of the order parameter. They have demonstrated 
that it is the magnetism which splits the critical temperature. Above the critical pressure, 
the hexagonal symmetry is restored. Experiments above this pressure have shown that there 
are still two phase transitions as a function of field - this means that the field direction itself 
couples to the internal degrees of freedom in the two-component order parameter. This 
only occurs in the E\ g picture, which appears to be the only choice fully consistent with all 
experiments. 

We would like to acknowledge useful discussions and correspondence with M. Norman, A. 
Garg, and particularly D. Cox. This work was supported by the National Science Foundation 
through grant no. DMR-9214739. 
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FIGURES 



FIG. 1. Consequences of the effective field in the s-wave case, (a) Eigenvalue curves for the 
original lattice with k E L. The lines are the solutions to Eq. 2. 14] . However, all transitions are 
suppressed by the effective field except the original H C 2 line. In the case of arrows, the repulsion 
of the boundary to unphysical values of H and T takes place. In the case of cross-hatching, the 
transition is converted to a crossover. The numbers in parentheses denote the class of the line, (b) 
Curves for the new lattice with k £ L' . As in (a), except that the new lattice interpenetrates the 
old one. All lines are repelled by the effective field, and the C n fc corresponding to these boundaries 
never become nonzero. 



FIG. 2. Consequences of the effective field in the d-wave case, (a) Eigenvalue curves for the 
original lattice with k € L setting the eigenvalues of the quadratic form in Eq. 2.28 equal to zero, 
neglecting the fourth-order terms. Most transitions are converted to crossovers by the effective 
field except the original H C 2 line. The numbers in parentheses denote the class of the line, (b) 
New lattice with k £ V . As in (a), but the new lattice interpenetrates the old one. Transitions 
corresponding to Co a and Co b are repelled only a short distance (single arrows). The dashed lines 
show the final positions of these boundaries after taking into account the effective field. A single 
internal transition line remains. 



FIG. 3. The flux lattice for rj y . The filled circles indicate the singularities in the flux lattice. 
The Wigner-Seitz primitive unit cell is indicated by the dashed line. Its basis vectors are also 
shown. The integral I2 (| I2 | = |< ( r lx) 2 'n'y >l) i s non-zero only if the singularities in the r\ x flux 
lattice are located at the b (open square), c (open oval) or d point (x), or one of the symmetrically 
equivalent points. These points for the b (c) point are the filled squares (filled oval). The axes are 
z = z/£ (|| c-axis) and y = y/£. The diagram is drawn to scale using the values of the parameters 
we obtain through our fits. The flux lattice for r/ x is identical except for greater elongation in the 
z direction. 



FIG. 4. Phase diagram when the field is in the basal plane. The data points are from ultrasonic 
velocity measurements and are taken from Ref. ||, Fig. 3. The solid lines are our fit without the 
renormalization discussed in the text. The dashed line coresponds to a renormalization of Q x by 
2.67 or Q x = 2.67Q X . The constants used to make these graphs are T x = 0.458 K, T y = 0.504 K, 
S C 2x = -9.26 T/K, S c2y = -4.39 T/K, a x /a = 0.0138 K/T 2 , and a d /a Q = 0.0193 K/T 2 



FIG. 5. Phase diagram when the field is in the z-direction. The lines are the theoretical fits to 
H C 2 (solid line) and the inner transition (dashed line). The data points are from ultrasonic velocity 
measurements and are taken from Ref. 0, Fig. 3. In addition to the constants used to make Fig. |] 
we have used a /K x = 5.6 x 10 12 K^cm" 2 , e = 5.26 x 10~ 5 G _1 , a z /a = 6.3 x lO" 10 K/G 2 , and 
(K 2 -K 3 )/K x = 0.1. We use {K 2 +K 3 )/K 1 = (S c2x /S c2y ) 2 -l + hca eAT[(S C 2 X /S C 2 y ) 2 + l)}/{2eK 1 ) 
and if 4 /iTi = [(1 + a eAT '/ 'Kx)[2eKiS c2x / '{hca^)} 2 ]- 1 to obtain (K 2 + # 3 )/ifi = 1.0 and 
K±/K x = 7.20. 
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FIG. 6. Susceptibility of UPt3 for fields oriented along the crystal's a-axis (circles), b-axis 
(triangles), and c-axis (squares). The graph is taken from Ref. Fig. 2.1. 



FIG. 7. Fermi surface of UPt3 from Ref. |36|, Fig. 5. The Fermi surface shown was calculated 
in the local density approximation using the Dirac-relativistic linear muffin-tin orbital method. 
The stripes show the j z content of the Fermi surface: dotted for \ j z \ = 1/2, right-hatched (/ / /) for 
\jz\ = 3/2, and left-hatched for \j z \ = 5/2. The \j z \ = 7/2 component is small over the entire Fermi 
surface. Note that the parts of the Fermi surface around the T and K points are predominantly 
\jz\ = 3/2 or 5/2, while the parts around the A point have well-mixed |j z |'s. Fig. 2(b) of Ref. |35| 
is similar. 



FIG. 8. The crossing of the H C 2 line when the field is the basal plane (solid line) and the H C 2 
line when the field is in the z-direction (dashed line). The data points are from ultrasonic velocity 
measurements and are taken from Ref. ||, Fig. 3. The diamonds are for the case when the field in 
the basal plane (H || ab) and the crosses are for the case when the field is in the z-direction (H \\ 
c). 



FIG. 9. Pressure dependence of the phase diagram with the field in the basal plane in the E\ g 
model. The phase diagram is plotted at pressures (P) of (a) P = 0, (b) P = Pn/2 (1-85 kbar), 
(c) P = Pjv (3.7 kbar), and (d) P = (3/2) (5.55 kbar). Here P/v is the pressure above which 
the temperature splitting vanishes. As discussed in the text the theoretical inner transition line 
for temperatures below the tetracritical point (H*(T)) has been renormalized. The data points in 
(a) are taken from Ref. ||, Fig. 3. The variation of the transition temperatures with pressure is 
taken from Ref. |41]. 



FIG. 10. Pressure dependence of the phase diagram with the field in the basal plane with 
(Ki + K^)/K\ = (E2u model). The phase diagram is plotted at pressures (P) of (a) P = 0, (b) 
P = (1/2)P N (1.87 kbar), (c) P = P N (3.7 kbar), and (d) P = (3/2)P N (5.55 kbar). In order to 
obtain a better fit in this model we have changed the values of some of our input parameters. In 
these phase diagrams a x = ad = 0, S C 2 X = —6.66 T/K, S C 2 y = —4.07 T/K, T x = 0.465 K, and 
T y = 0.509 K. The data, renormalization of H*(T), and all other input parameters are the same 
as those for Fig. ^. 
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K. A. Park and Robert Joynt: Fig 1(a) 




K. A. Park and Robert Joynt: Fig 1(b) 




K. A. Park and Robert Joynt: Fig. 2(b) 




K. A. Park and Robert Joynt: Fig. 3 
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K. A. Park and Robert Joynt: Fig. 7 
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K. A. Park and Robert Joynt: Fig. 9(a) 
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K. A. Park and Robert Joynt: Fig. 10(a) 
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K. A. Park and Robert Joynt: Fig. 10(c) 



